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ABSTRACT 

We study the chemical properties of the stellar populations in eight simulations of 
the formation of Milky- Way mass galaxies in a ACDM Universe. Our simulations 
include metal-dependent cooling and an explicitly multiphase treatment of the effects 
on the gas of cooling, enrichment and supernova feedback. We search for correlations 
between formation history and chemical abundance patterns. Differing contributions 
to spheroids and discs from in situ star formation and from accreted populations 
are reflected in differing chemical properties. Discs have younger stellar populations, 
with most stars forming in situ and with low a-enhancement from gas which never 
participated in a galactic outflow. Up to 15 per cent of disc stars can come from 
accreted satellites. These tend to be a-enhanced, older and to have larger velocity 
dispersions than the in situ population. Inner spheroids have old, metal-rich and a- 
enhanced stars which formed primarily in situ, more than 40 per cent from material 
recycled through earlier galactic winds. Few accreted stars are found in the inner 
spheroid unless a major merger occurred recently. Such stars are older, more metal- 
poor and more a-enhanced than the in situ population. Stellar haloes tend to have low 
metallicity and high a-enhancement. The outer haloes are made primarily of accreted 
stars. Their mean metallicity and a-enhancement reflect the masses of the disrupted 
satellites where they formed: more massive satellites typically have higher [Fc/H] and 
lower [a/Fe]. Surviving satellites have distinctive chemical patterns which reflect their 
extended, bursty star formation histories. These produce lower a-enhancement at given 
metallicity than in the main galaxy, in agreement with observed trends in the Milky 
Way. 

Key words: galaxies: formation, galaxies: evolution, galaxy: stellar content, cosmol- 
ogy: theory 



1 INTRODUCTION 

According to the current cosmological paradigm, galaxy for- 
mation occurred through the condensation of gas in the cores 
of a hierarchically growing population of dark haloes. This 
strongly nonlinear process can be followed in detail only with 
sophisticated numerical simulations. Advances in numerical 
techniques and in computational power, together with the 
accumulation of more precise and more detailed observa- 
tions of the galaxy population, have enabled substantional 
improvements in our understanding of cosmic structure for- 
mation. Nevertheless important open questions remain, in 
particular on galactic scales ( Mo et al.|20 10). Many of these 
issues are related to the evolution of the baryons. A vari- 
ety of physical processes, acting on different scales and with 
differing efficiency, are believed to play a role in determin- 



ing the properties of galaxies. Direct simulation of galaxy 
formation in a ACDM context has evolved greatly in recent 



years (e.g. 


Springel 2005 


Scannapieco et al. |2006 Gover- 


|nato et al. 


2007; Brook et al. 2007 Stinson et al. 2006 


Sijac 


ri et a 


. 2007 Scannapieco et al. 


2009| |Agertz et al. 


[2009 
[2011 


Schaye et al.| 2010 Agertz et al. 


|2011| |Brook et al. 


with improved descriptions of cooling, star formation, 



black hole formation, and luminous, mechanical and mate- 
rial feedback from star-forming regions and active galactic 
nuclei. Nevertheless, most aspects of this critical baryonic 
physics are still described by schematic, phenomenological 
sub-grid models, since it is not yet possible to achieve the 
small-scale resolution needed to simulate them ab initio. 

Confronting such simulations with observation is a key 
route to improving the models and developing a better 
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understanding of galaxy growth in its cosmological con- 
text. Chemical patterns can be powerful tools for advanc- 
ing this programme. The elements are synthesized in stars 
and ejected into the interstellar medium (ISM) through stel- 
lar winds and supernovae (SNe). Different SN types release 
their nucleosynthesis products on different timescales, im- 
printing characteristic chemical patterns on the stellar pop- 
ulations which form from the gas they enrich (as observed 
in the Milky Way, see e.g|Bland-Hawth orn fc Freeman|2003"l 
WysepoTo] ). SN products are mixed with other interstellar 
gas by internal mechanisms, or may even be thrown out of 
the galaxy as part of a wind. Environmental effects such as 
interactions with other galaxies or ram pressure stripping as 
a result of motion through the intergalactic medium, also 
play a role by mixing new heavy elements into the inter- 
galactic medium. Galaxy mergers induce major starbursts 
which substantially modify the abundance patterns in the 



mass haloes from the Aquarius Project ( Springel et al. 2008 1 



ISM (e.g. Perez et al. 2006 Pi Matteo et al. 2009 Rupke 



et al. 2010 Perez et al. submitted). Even minor mergers 
with infalling satellites modify metal abundance patterns 
by contributing stars and gas with distinctive abundances 



and by triggering new star formation activi ty (e.g. |Barnes 
fc Hernquist|1996||Mihos fc Hernquist|1996| |Tissera|2000[ ). 

The possible contribution of disrupted satellites to 



the stellar populations of thin discs (Abadi et al. 20031, 
thick discs (Huang & Carlberg 1997[ ) and bulges ( |Rahimi| 
et al. 2010 1 has been much debated. Some simulations 



have suggested that stellar haloes may have formed primar- 



ily through satellite disruption (e.g. Johnston et al. 2008 



Cooper et al. 20101, but there have also been claims for 



a dual formation route, depending on the early merging 



history of each individual galaxy (e.g. Zolotov et al. 2009 



20101; this may be supported by abundance patterns in the 



stellar halo and satellites of the Milky Way (Venn et al 



20041. Different evolutionary paths are expected to leave 
their fingerprints on the chemical and dynamical properties 
of the stellar populations, and analysing such fingerprints 
should help us to learn more about how galaxies formed. 

Much work has been devoted to understanding if chem- 
ical patterns can indeed be related to specific events during 
galaxy formation (e.g. |Molla fc Ferrini|19 95 Chiappin i et al.| 
19971. Within the current hierarchical paradigm, the nonlin- 



earity and complexity of structure growth require chemical 
and structural evolution to be treated simultaneously and 
consistently if reliable conclusions are to be drawn. In the 
last decade a variety of treatments of chemical evolution 
have been implemented in the cosmological hydrodynam- 
ics codes used to simulate galaxy formation (e.g. |Raiteri| 



eTaT1^996l |Mosconi et"aLl[200Tj |Lia et al.|[2002 


Kawata &| 


Gibson 2003; Kobayashi et al.||2007| |Wiersma et al.| 2009 


Martinez-Serrano et al. 


2008 


1. These typically adopt the 



standard metal production assumptions from analytic mod- 
els of chemical evolution (e.g. |Tinsley fc Larson] [!979}|MaFl 
teucci fc Greggio||1986[) but follow the formation and trans- 



port of elements among the various components of the sys- 
tem in a dynamically consistent way. Such codes have al- 
ready given a number of encouraging results, reproducing 
general trends observed in the galaxy population such as 



the mass-metallicity relation, for example (Governato et al. 
2007| |Kobayashi et al.|2007l ) 



The simulations analysed in the current paper follow 
the formation of galaxies within a set of eight Milky- Way- 



The resulting central galaxies have been extensively studied 



in papers by|Scannapieco et al.| (|2009), Tisse ra et al.| (2010), 

(20TT|. 



Scannapieco et al. 



Scannapieco et al. 



( 2010| ) and 

These studies highlighted the difficulty in forming spirals 
similar to the Milky Way within the ACDM cosmology (but 
see Guedes et al. (2011) who were able to produce an L* 
disk galaxy with a low i-band bulge-to-disk ratio in a single 
high-resolution simulation which assumed a higher surface 
density threshold for star formation than earlier simulations 



like our own.). Scannapieco et al. ( 2009 1 found that although 
most of the simulated galaxies had a substantial disc at some 
stage during their evolution, these discs are often disrupted 
by interactions with satellites or destabilized by misaligned 
inflowing cold gas. These processes feed the growth of stel- 
lar spheroids, which end up as the dominant components by 
mass in all the simulations. Thus, at z = none of our sim- 
ulated haloes hosts a spiral with a disc-to-total mass ratio 
similar to that of the Milky Way. In spite of this, the final 
discs and spheroids show general dynamical properties which 
are quite consistent with those of observed spiral galaxies. 
While the spheroids are dominated by old, dynamically hot 
stars (and tend to rotate more slowly than observed bulges), 
the discs have the youngest stars, they formed typically from 
the inside out, and they are dynamically cold. Residual cold 
gas forms a well-behaved disc in equilibrium within its po- 



tential well ( |Scannapieco et ah 20111. The fact that these 
two components show distinct star formation histories and 
distinct dynamical properties prompted us to explore their 
chemical properties, seeking chemical signatures which can 
be related to specific details of their formation histories. 

The simulations studied here use a specially enhanced 
version of GADGET-3 which is itself an optimisation of 



GADGET-2 (Springel 20051. Our version includes a chem- 
ical and energy feedback model which is consistently cou- 
pled to a unique multiphase representation of the interstel- 
lar medium (ISM) which allows collocation, material ex- 
change, and relative motion between the different phases. 
This scheme allows us to follow the dynamical and chemi- 
cal properties of baryons throughout the galaxy formation 
process. Dense and cold gas is transformed into stars which 
can explode as SN. These inject entropy into some of the 
remaining cold gas, converting it to the hot diffuse phase 
and driving powerful, mass-loaded winds. These thermody- 
namic changes are generated on a particle-by-particle basis 
in a way which is consistent with local astrophysical condi- 
tions. As a result star-forming regions drive winds without 
the need to temporarily suspend cooling, to impart large 
kicks to particles, or to suppress hydrodynamical interac- 
tions with surrounding particles. This naturally results in 
self-regulated winds which are adapted to the depth of the 
potential well in which star formation is occurring (jScanna- 
pieco et al.||2006l [2009] |Sawala et aLl|20Tol [20TT| |E>e Rossi 



et al.|20T0 l. In this framework, gas is progressively enriched 



as subsequent stellar generations release new elements into 
the ISM, so that the chemical structure of each galaxy is 
intimately related to the formation paths of its various com- 
ponents. In this paper, we focus on this relation between 
chemical properties and assembly history. In a subsequent 
paper, we will analyse how gas flows lead to abundance gra- 
dients (Tissera et al. 2011 in preparation). 

In Section 2 we describe our simulations and the main 
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aspects of our chemical model. We also include a subsection 
explaining the criteria we adopt to separate our galaxies 
into their dynamical components. In Section 3 we study the 
chemical properties of these components, together with their 
formation histories. Section 4 studies the chemical proper- 
ties of stars formed in situ and contrasts them with those 
of stars accreted from tidally disrupted satellites. Section 5 
compares the mean chemical abundances of stars from dis- 
rupted satellites with those of stars in surviving satellites. 
We summarize our findings in the concluding section. 



2 THE SIMULATIONS 

We have analysed galaxy formation within a set of eight 
haloes selected and simulated at high resolution as part 
of the Aquarius Project (Springcl ct al. 2008[). These ob- 
jects were identified in a lower-resolution version of the 



Millennium-II Simulation ( Boylan-Kolchin et al. 20091, a 
cosmological simulation of a cubic region 100 Mpc h _i on 
a side, imposing a mild isolation criterion and requiring a 
mass similar to that of the Milky Way's halo at z = 0. In 
the original project these haloes were simulated at a variety 
of resolutions following the dark matter only. For the simu- 
lations discussed here, gas particles were added to the initial 
conditions on a cubic mesh displaced by half the interparti- 
cle separation from the dark matter particle mesh. Details 



can be found in Scannapieco et al. ( 2009 1 



All these simulations assumed a ACDM cosmogony with 
n m = 0.25, Q A = 0.75, n h = 0.04, cr 8 = 0.9, n s = 1 and 
H = 100 h kms^Mpc^ 1 with h = 0.73. Our galaxy for- 
mation simulations were carried out with maximum gravi- 
tational softenings in the range eg = 0.5 — 1 kpc h~ , dark 
matter particle masses of order 1O 6 M0/i _1 and initial gas 
particle masses of order 2 x 1O 5 M0/i _1 . Since their virial 
masses are in the range 5 to 11 x 10 1:l Mq/i _:l , the final 
systems all have about a million particles in total within 
the virial radiuiQ In Table [2] we summarize their princi- 
pal characteristics. Our naming convention (Aq-A-5, etc.) 
follows that of earlier papers. The capital letter (A to H) 
designates the particular halo considered, whereas the num- 
ber (1 to 5) indicates the resolution level of the simulation 
according to the original Aquarius convention. All simula- 
tions studied in this paper are at level 5 the lowest of the 
original Aquarius resolutions. Further details can be found 
kill 



2.1 The galaxy formation simulation code 

The version of GADGET-3 used to carry out these simula- 
tions has been extended to include an explicit multiphase 
model for the gas component which allows a cold, dense, 
possibly star-forming phase and a hot, diffuse phase to in- 
terpenetrate, to move relative to each other, and to ex- 
change material. This implementation also includes metal- 
dependent cooling, star formation, and chemical and ther- 
mal feedback into both gas phases, as described in |Scan-| 
|napieco et al.| ( |2005[ ) and [Scannapiec o et ak] ( |2006| ) . This 
scheme has been successful in reproducing the star forma- 
tion activity of galaxies during both quiescent and starburst 
phases, and can drive substantial mass-loaded galactic winds 
from star-forming galaxies with a wide range of masses at 
speeds that consistently reflect the depth of the individual 
potential wells (Scannapi eco et aTj 2006, 2008; Sa wala et al.| 



2010| |2011| |De Rossi et al.||2010[ ). The scheme does not re- 
quire any ad hoc "fixes" (for example, temporary suppres- 
sion of radiative cooling and/or of hydrodynamic coupling 
for wind particles, or the imposition of discontinuous mo- 
mentum changes to "kick-start" a wind) to ensure efficient 
wind generation in SPH simulations, and there are no scale- 
dependent parameters. As a result, it is particularly well 
suited for the study of galaxy formation in a cosmological 
context where galaxies with a wide range of properties form 
simultaneously. The version used here includes chemical en- 
richment by Type II and Type la supernovae (SNII and 
SNIa, respectively) as we describe next. 



2.1.1 The Chemical Model 

The chemical evolution model adopted in this work is that 



developed by |Mosconi et al. (20011 and adapted by Scan- 
|napieco et al.| (|2009|> to GADGET-3. We consider chemical 



elements which are synthesised in stellar interiors and then 
ejected into the ISM by SNII and SNIa events. A Salpeter 
Initial Mass Function is assumed with lower and upper mass 
cut-offs at 0.1 M Q and 40 M , respectively. We follow 12 dif- 
ferent chemical isotopes: H, 4 He, 12 C, ie O, 24 Mg, 28 Si, 56 Fe, 
14 N, 20 Ne, 32 S, 40 Ca and 62 Zn. Initially, all baryons are in 
the form of gas with primordial abundances, i.e. 76% H and 
24% 4 He by mass. 

SNII are considered to originate from stars more mas- 
sive than 8 Mq, adopting the metal-dependent yields of 



Scannapieco et al. (2009) 



Woosley & Weaver (1995). We use the metal- and- mass de- 



Scannapieco et al j (120091) andlScannapieco et all (120111) pendent lifetimes given by the fitting formulae of jRaitcri 



studied the formation of the disc and the spheroidal compo- 
nents of these systems and their dynamical and photometric 



et al. (19961. For SNIa, we adopt the W7 model of Thiele- 



mann et al. ( 1993 1 which assumes that SNIa events originate 



properties, while Tissera et al. (20101 analysed the reponse 



of their dark matter haloes to the galaxy formation pro- 
cess. The main galaxies have a variety of morphologies with 
varying bulge-to-disc ratios and strong bars in several cases. 
Nevertheless, all of them are bulge dominated, and one (Aq- 
F-5) has no detectable disc at all at z = 0. In this paper, we 
focus on the chemical properties and the assembly histories 
of the various stellar components of the galaxies. 



1 We define the virial radius (^200) as the radius within which 
the mean mass density of all components is 200 times the critical 
value. The virial mass is then the total mass within this radius. 



from CO white dwarf systems in which mass is transferred 
from the secondary to the primary star until the Chan- 
drasekhar mass is exceeded and an explosion is triggered. 
For simplicity, we assume that the lifetime of the progen- 
itor systems is distributed uniformly and at random over 
the range [0.7,2] Gyr. To calculate the number of SNIa, we 
adopt an observationally motivated rate of 0.3 relative to 
SNII ( |Mosconi et aT^|200Tt . 

The ejection of chemical elements is grafted onto the SN 
feedback model so that chemical elements are distributed to 
the cold and the hot gas phases surrounding a given star 
particle. The fraction of elements going to each phase is 
regulated by a free parameter e c which also determines the 
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amount of SN energy received by each phase. While the in- 
jection of energy follows two different paths depending on 
the thermodynamical properties of the gas, chemical ele- 
ments are added to both phases at the time each SN oc- 
curs. All the simulations of this paper have been run with 
e c = 0.5. Some tests of the effect of varying this parameter 



can be found in Scannapieco et al. ( 2006 1 



2.2 



The stellar components of the simulated 
galaxies 



As shown by Scannapieco et al. (20091, the simulated 



galaxies have a variety of morphologies and contain both 
spheroidal and disc components. Our main contribution to 
understanding how these galaxies formed comes from linking 
their dynamical evolution with specific chemical patterns in 
their stellar populations. 

Several methods for separating spheroids and discs have 



been proposed (e.g. Abadi et al. 2003 Governato et al 



2009| |Scannapieco et al.||2009 l. In this work we use a pro 



cedure similar to that of Abadi et al. (20031. We measure 
e = J z j J Zt max{E) f° r each star, where J z is the angular 
momentum component perpendicular to the disc plane and 
Jz,max(S) is the maximum J z over all particles of given to- 
tal energy, E. A star on a prograde circular orbit in the 
disc plane has e ~ 1. In Fig. [T] we show e as a function of 
total energy for the stars in each of the Aquarius haloes. 
The colour code indicates stellar metallicity, [Fe/H]. Discs 
are easily identified as star concentrations near e = 1 while 
stars in the spheroidal components usually scatter around 
e = 0. The bulge is located at the lowest (most bound) en- 
ergies, while the outer halo is at the highest (least bound) 
energies. 

We define stars to belong to different dynamical com- 
ponents according to their location in Fig[l] We consider 
stars with e > 0.65 to be part of a disc. This value is small 
enough to include stars in both thin and thick discs. Par- 
ticles which do not satisfy this requirement are taken to 
be part of the spheroid. Motivated by observations of the 
Milky Way spheroid which report differences in stellar kine- 
matics and chemical abundances as one moves outwards 



(Carollo et al. 2007 Zoccali et al. 20081, we separate our 



spheroids into three subregions according to their binding 
energy. The central spheroid ("bulge") is defined to consist 
of stars more bound than the minimum energy (£J ce n) of 
stars with r > 0.5 x r op t where the characteristic radius 
r opt is defined to enclose 83% of the baryonic mass of the 
simulated galaxy. Central spheroid stars thus all have or- 
bital apocentres less than about 0.5r opt . Stars more weakly 
bound than -Been but with E < Ei nnol , the minimum energy 
of all stars with r > 2 x r op t are considered to form the 
inner halo. Stars with E > -Binncr are defined to belong to 
the outer halo and thus all have orbital pericentres outside 
about 2r op t- Although there is some arbitrariness in these 
criteria, the main features of the components change only 
slowly with the reference energies. Our criteria are chosen 
so that the definition of the spheroidal components adapts 
to the overall size of each individual galaxy. Fig. [I] indicates 
the two bounding energies (E ccn and i5i nncr ) for each Aquar- 
ius system. 

The disks extend over a relatively small range of ener- 
gies and thus radii in Aq-A and B, but over a more sub- 



stantial range in Aq-C, D, E and G. Aq-F and H appear to 
have very little disc at all. The relative mass of the discs is 
seen more easily in Fig[2] where we give histograms of e for 
all stars within 2r opt in each of our galaxies. Returning to 
Fi^T] the central spheroid is seen to have rather little rota- 
tion in Aq-A, C, D, F, G and H, but is clearly rotating in the 
same sense as the disc in Aq-E. In Aq-B the inner spheroid 
is rotating in the same sense as the disc in the innermost 
regions but in the opposite sense at intermediate energies. 

The inner stellar halo is essentially non-rotating in all 
cases except Aq-E and F. In the former there is no clear 
separation of the inner halo and disk populations, as is par- 
ticularly evident in Fig[2] The outer halo is non-rotating in 
all cases except Aq-F, and in most galaxiees it shows stellar 
concentrations associated with satellite galaxies and stellar 
streams (i.e. disrupted satellites). For our main analysis we 
remove all satellites detected as individual bound structures 
in order to focus on the "diffuse" stellar halo. The chemical 
abundances of surviving satellites are analysed separately 
in the last section of this paper. Colour variations between 
components in different panels of Fig[l]indicate different lev- 
els of chemical enrichment. In general, strongly bound stellar 
populations are more chemically enriched than those found 
at lower energies. However, disc stars are clearly more chem- 
ically enriched than spheroid stars of similar energy. 



3 FORMATION HISTORIES OF THE 

STELLAR COMPONENTS AND THEIR 
GLOBAL CHEMICAL PROPERTIES 

It has long been suggested that particular chemical patterns 
are linked to key events in the formation history of galaxies. 
In this section, we will address this question by comparing 
formation history with 2 = abundance pattern for the 
various components of our galaxies. 

In Fig. [3j we present star formation histories for indi- 
vidual stellar components in the form of cumulative plots 
of the fraction of the final stellar mass formed as a func- 
tion of age. In such plots, a step reflects a starburst. The 
central spheroid and the inner and outer haloes have quite 
similar SF histories, with most of their stars formed at very 
high redshift. The oldest stars are split between the central 
spheroid and the haloes. The discs contain younger stellar 
populations which form in a series of episodes spread over 
the Hubble time. In the case of Aq-F-5, no disc component 
survives at z = 0. 

Star formation histories vary from system to system as a 
consequence of differing evolutionary paths, leading to vari- 
ations in chemical properties between our galaxies and be- 
tween components within a galaxy. This can be seen in Fig. 
[4]where we show [Fe/H] distributions for the individual stel- 
lar components in four of our simulated haloes. Interestingly, 
while the central spheroids, the inner and the outer stellar 
haloes all have similar star formation histories (see Fig. [3| 
their [Fe/H] distributions are very different. Comparing in- 
dividual stellar components between galaxies, we also note 
significant variations in the abundance distributions. 

The relative abundance of a-elements with respect to 
Fe is often identified as a key property to characterise the 
formation histories of stellar populations. We have measured 
[O/Fe] distributions for the individual stellar components in 
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Table 1. General characteristics of our main galaxies and their haloes at z = 0. The left column shows the halo name. 7*200 and M200 are 
the virial radius and mass, respectively. N^m, N gBls , N^ ot and AT* are the total numbers of dark matter, gas, hot gas and star particles 
within T200 1 respectively. r op t is the optical radius of the galaxy and M° pt and N° pt are the stellar mass and the number of star particles 
estimated within twice r opt . 



Halo 


»*200 
h~ x kpc 


M200 
lO 12 /i- 1 M 


Mhot 

io 10 /»- 1 m s 


M co l d 

lO lo ft -1 M0 


AW 


AT gas 


N hot 


N* 


r op t 
h kpc 


M opt 

10 10 /» -1 M a 


N° pt 


Aq-A-5 


169.42 


1.10 


4.1 


1.1 


529110 


92159 


72157 


333409 


13.1 


5.92 


277382 


Aq-B-5 


132.10 


0.52 


1.2 


0.6 


435330 


63580 


41312 


291313 


13.0 


2.53 


244511 


Aq-C-5 


173.19 


1.18 


3.4 


0.6 


681143 


97796 


83845 


548964 


11.7 


5.93 


432232 


Aq-D-5 


170.63 


1.09 


1.7 


0.2 


599438 


82736 


75741 


377866 


10.8 


4.41 


287734 


Aq-E-5 


149.93 


0.79 


2.1 


0.6 


554245 


80524 


63458 


512809 


10.1 


4.97 


406399 


Aq-F-5 


142.74 


0.67 


1.6 


0.2 


680129 


79345 


71124 


699052 


10.3 


5.39 


525101 


Aq-G-5 


142.64 


0.68 


1.2 


0.5 


679177 


74675 


52828 


387751 


10.3 


5.63 


334108 


Aq-H-5 


132.94 


0.53 


0.5 


0.1 


515392 


23133 


20671 


546482 


7.60 


4.73 


425680 



Table 2. Properties of the stellar components of the main galaxies at z = 0. M^ ac , M ccn , Mi nncr and Aioutcr are the stellar masses 
of disc, central spheroid, inner halo and outer halo of the main Aquarius galaxies as obtained from the kinematic decomposition. All 
masses are given in units of 10 9 /i — *Mq. D/T 1 and D/T 2 are the disc-to-total stellar mass considering as total mass only the stellar discs 
and the central spheroids (T 1 ) and also the stellar inner haloes (T 2 ). The last column show D/TP hoto photometric ratios calculated by 
|Scannapieco et al.|{20l6") using an observationally based approach (for Aq-F-5 amd Aq-H-5, they could not obtain a reliable photometric 
decomposition) . 



Halo 


M diBC 


M con 


-^Inncr 


A^Outer 


D/T 1 


D/T 2 




Aq-A-5 


3.2 


3.4 


14.6 


2.4 


0.08 


0.06 


0.32 


Aq-B-5 


2.2 


16.9 


6.4 


1.4 


0.11 


0.09 


0.42 


Aq-C-5 


1.6 


37.8 


14.2 


3.6 


0.29 


0.23 


0.49 


Aq-D-5 


8.9 


24.0 


13.6 


4.0 


0.27 


0.19 


0.68 


Aq-E-5 


1.1 


25.2 


15.3 


2.1 


0.30 


0.21 


0.40 


Aq-F-5 




39.0 


9.5 


0.5 








Aq-G-5 


6.0 


16.0 


7.3 


0.7 


0.27 


0.20 


0.60 


Aq-H-5 


1.3 


15.6 


22.3 


2.5 


0.07 


0.03 





four of our galaxies as displayed in Fig. [4] (upper raw). 
As in the Milky Way, discs show the lowest level of a- 
enhancement. This is because their stars are younger than 
the other simulated stellar populations and formed over an 
extended period of time. As a result, the gas from which they 
formed had time to be enriched by SNIa. Stars in the central 
spheroids are a-enhanced but have the highest [Fe/H] values 
of all the populations, reflecting the fact that they formed 
in a short and strong starburst and were able to retain a 
large fraction of their SN ejecta. These general properties 



are in agreement with data for the Milky-Way bulge (Zoc- 
cali et al.|20 08). The largest a-enhancements in our models 
are, however, found in the inner and outer stellar haloes. 

Although the typical [Fe/H] values found for our bulge 
and disc components are lower than for their Milky Way 
counterparts, the trends between components are simi- 
lar to those observed. There are clear variations between 
the various Aquarius galaxies, both in [Fe/H] and in a- 
enhancement, but the above mentioned trends seem to be 
common to all of them. Remaining differences with the 
Milky Way (for example the G-dwarf problem evident in 
Fig. [4| must, at least in part, be due to the fact that 
all the Aquarius galaxies have disc-to-spheroid mass ratios 
much smaller than that of the Milky Way. We also note that 
the metallicity distributions have large standard deviations 
which might suggest the need for more efficient metal mix- 



ing {Martinez-Serrano et al.|2008||Wiersma et al.|2009[ ). For 
the discs and inner haloes, these large dispersions make it 
difficult to distinguish chemical differences. Nevertheless, for 
the best discs, Aq-C-5, Aq-D-5 and Aq-G-5, a Kolmogorov- 
Smirnov test shows the metallicity distributions of the discs 
and inner haloes to differ at confidence levels between 70 
and 94 percent. Only for Aq-B-5 and Aq-H-5 do the the 
metallicity distributions of these components appear indis- 
tinguishable. In general, for all our galaxies, the spheroidal 
components show a decrease of [Fe/H] as one moves out- 
wards. Their stellar populations are more a-enriched and 
older than those populating the corresponding disc compo- 
nents. Stars in the haloes are, on average, more a-enriched 
than those in the central spheroids. 

To test how numerical resolution affects our estimates 
of chemical abundances, we have analysed two lower resolu- 
tion versions of Aq-E which we label Aq-E-6 and Aq-E-7^] 
Within r2oo these have totals of about 160,000 and 80,000 
particles, respectively, whereas Aq-E-5 has more than one 
million. These two lower resolution simulations were also 



analysed by. Scannapieco et al. (2009 20111 and Tissera 



et al. (20101 who reported good agreement between their 



dynamical and structural properties and those of Aq-E-5, 



Simulations Aq-E-6 and Aq-E-7 correspond to those labeled as 



Aq-E-6b and Aq-E-6 in Scannapieco et al. (2011) 
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Table 3. Chemical properties of the stellar discs. The first column gives the encoding name of the halo. F pr0 g and F prom are the fractions 
of stars formed in the progenitor systems and of stars which have been part of a galactic wind. a/V is the ratio of velocity dispersion 
to mean circular velocity. Age is the mean stellar age in Gyr. [Fo/H] and [O/Fe] are median logarithmic abundances relative to solar. 
Standard deviations are included within parentheses. The superscripts "prog" and "sat" denote values corresponding to stars formed in 
situ and to stars acquired from disrupted satellites, respectively. 



Halo 


■fprog 


■fprom 


cr/yprog 


(cr/V) sat 


[Fe/H]P r °s 


[Fe/H] sat 


[0/Fe]P r °s 


[0/Fe] sat 


Age prog 


Age sat 


Aq-A-5 


0.97 


0.08 


0.86 


1.13 


-1.10(0.43) 


-1.75(1.35) 


-0.163(0.133) 


0.343(0.052) 


5.89(7.53) 


12.66(0.13) 


Aq-B-5 


0.91 


0.09 


0.67 


0.89 


-1.25 (0.62) 


-1.55(1.03) 


-0.174(0.154) 


0.094(0.288) 


4.08(4.19) 


10.63(4.92) 


Aq-C-5 


0.99 


0.09 


0.59 


0.76 


-0.86 (0.33) 


-1.37(0.81) 


-0.131(0.121) 


0.300(0.113) 


8.69(5.07) 


12.39(3.07) 


Aq-D-5 


0.89 


0.15 


0.68 


0.77 


-0.85 (0.34) 


-0.97(0.73) 


-0.096(0.114) 


0.192(0.142) 


7.99(6.47) 


10.81(1.26) 


Aq-E-5 


0.62 


0.30 


0.54 


0.60 


-0.95 (0.36) 


-0.83(0.58) 


-0.049(0.123) 


0.031(0.154) 


8.53(2.98) 


11.09(0.79) 


Aq-F-5 
Aq-G-5 


0.97 


0.12 


0.62 


0.73 


-0.94(0.38) 


-1.75(1.21) 


-0.092(0.109) 


0.076(0.259) 


8.06(7.07) 


12.22(2.85) 


Aq-H-5 


0.84 


0.21 


1.08 


1.02 


-0.90(0.39) 


-1.25(1.08) 


-0.137(0.131) 


0.241(0.151) 


7.72(7.03) 


12.36(0.59) 



Table 4. Chemical properties of the central spheroids. The first column gives the encoding name of the halo. -Fprog and F prom are the 
fractions of stars formed in the progenitor systems and of stars which have been part of a galactic wind, a is the mean velocity dispersion 
(in km s — x ), while Age is the mean stellar age (in Gyr). [Fe/H] and [O/Fe] are median logarithmic abundances relative to solar. Standard 
deviations are included within parentheses. The superscripts "prog" and "sat" label values corresponding to stars formed in situ and to 
stars acquired from disrupted satellites, respectively. 



Halo 


Fprog 


-Pprom 


[Fe/H]P r °s 


[Fe/H] sat 


[O/Fc] 


[0/Fe] sat 


Ag e P r °s 


Age sat 


Aq-A-5 


0.96 


0.43 


-0.61(0.25) 


-1.20(0.83) 


0.023(0.126) 


0.358(0.016) 


11.46(1.67) 


12.63(0.09) 


Aq-B-5 


0.79 


0.64 


-0.66(0.36) 


-0.92(0.55) 


-0.025(0.111) 


0.159(0.092) 


9.43(1.52) 


10.93(1.56) 


Aq-C-5 


0.95 


0.44 


-0.61(0.34) 


-1.12(0.56) 


0.038(0.122) 


0.353(0.033) 


11.82(1.18) 


12.52(0.04) 


Aq-D-5 


0.90 


0.60 


-0.62(0.29) 


-0.99(0.63) 


-0.003(0.110) 


0.214(0.084) 


10.75(1.39) 


11.95(0.20) 


Aq-E-5 


0.71 


0.75 


-0.40(0.31) 


-0.57(0.41) 


-0.031(0.108) 


0.014(0.109) 


10.56(3.95) 


11.53(0.65) 


Aq-F-5 


0.55 


0.53 


-0.52(0.32) 


-0.65(0.50) 


-0.061(0.103) 


0.037(0.118) 


10.58(4.62) 


10.93(3.02) 


Aq-G-5 


0.74 


0.73 


-0.58(0.35) 


-0.70(0.50) 


0.048(0.112) 


0.122(0.085) 


11.37(2.00) 


11.08(1.53) 


Aq-H-5 


0.83 


0.78 


-0.40(0.33) 


-0.90(0.46) 


-0.026(0.122) 


0.358(0.028) 


11.62(1.51) 


12.45(0.09) 



both for baryons and for dark matter. We find similarly good 
convergence for the chemical properties of each of our stel- 
lar components, as can be appreciated from Fig. 5 where we 
compare [Fe/H] distributions in the three simulations. 



3.1 Assembly histories 

In order to understand if the stars in each component were 
formed in situ (i.e. in the main progenitor of the final galaxy) 
or in smaller objects which were subsequently accreted, we 
follow each star from its time of formation to z = 0. This 
allows us to estimate the fraction of stellar mass in each 
component which formed in situ. We consider all self-bound 
bound clumps of more than 20 simulation particles to be 
separate objects, so that any star which formed in a clump 
other than the main progenitor is considered to have been 
accreted. 

Cumulative in situ and accreted fractions are shown in 
Fig. [6]as a function of the redshift at which the stars formed. 
The four panels give results for our four stellar components, 
and within each panel solid lines give in situ fractions while 
dashed lines give accreted fractions. Different colours refer 
to the different Aquarius galaxies as noted in the caption. 
Most of the stellar mass in the central spheroids was formed 
in situ varying from ~ 60 to ~ 95 per cent. Only in Aq-F- 
5, which underwent a major merger at z ~ 0.6, are the in 
situ and accreted fractions comparable. In contrast, in the 



outer stellar haloes the majority of the stars are accreted in 
all eight galaxies. (Remember that we have excluded stars 
which remain part of a bound satellite at z = when compil- 
ing these statistics.) The inner haloes show an intermediate 
behaviour with both fractions ranging between 25 and 75 
per cent. About 75 per cent of the inner halo of Aq-F-5 is 
accreted, most of it as part of the second object in the major 
merger. The z = discs are heavily dominated by in situ 
stars, with a maximum accreted contribution of ~ 15 per 
cent. Note however, that the largest accreted fractions are 
for Aq-E-5, where Figs [I] and [2] show there to be no clear 
separation between disc and inner halo, and Aq-H-5, where 
the disc component is very weak. 

To correlate assembly histories with chemical abun- 
dances, we compare the fraction of stars formed in situ 
(fprog) in each stellar component to the median abundance 
of the stars at z = 0. In Fig. [7] we show medians of [Fe/H] 
as a function of fprog (upper panel). As can be seen, stellar 
components become more metal-rich as one moves inward; 
the central spheroid is always much more chemically en- 
riched than the corresponding outer halo. The trends are 
encouragingly consistent with Milky Way data. The median 
[Fe/H] of central spheroids appears independent of in situ 
fraction, whereas in the outer halo is there a clear trend for 
higher [Fe/H] to correspond to smaller in situ fraction. 

In the lower panel of Fig. [7] we plot median [O/Fe] 
against Fprog- The weak a-enhancement of disc stars is a 
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Table 5. Chemical properties of the inner diffuse stellar haloes. The first column shows the encoding name of the simulations. Fp r0 g 
and -Fp r om are the fractions of stars formed in the progenitor systems and of stars which have been part of a galactic wind, a and Age 
are the mean velocity dispersion (in km s _1 ) and mean age (in Gyr), respectively. [Fe/H] and [O/Fe] are median abundances. Standard 
deviations are included within parentheses. The labels "prog" and "sat" denote the values corresponding to stars formed in situ and to 
stars acquired from satellites. 



Halo 


p 

-fprog 


p 

L prom 


[Fe/H]P r °s 


[Fe/H] sat 


[0/Fe]P r °e 


[C-/Fe] sat 


AgcP r °s 


Age sat 


Aq-A-5 


0.75 


0.15 


-1.17(0.52) 


-1.45(1.07) 


0.060(0.145) 


0.351(0.054) 


10.40(4-05) 


12.60(0.15) 


Aq-B-5 


0.32 


0.27 


-1.37(0.85) 


-1.24(0.79) 


0.031(0.150) 


0.113(0.223) 


10.45(1.94) 


10.12(2.33) 


Aq-C-5 


0.56 


0.26 


-1.05(0.81) 


-1.20(0.8) 


0.002(0.141) 


0.311(0.008) 


11.50(0.92) 


12.37(0.40) 


Aq-D-5 


0.57 


0.27 


-1.09(0.54) 


-1.03(0.75) 


-0.079(0.135) 


0.152(0.146) 


9.96(2.67) 


11.47(0.53) 


Aq-E-5 


0.47 


0.43 


-0.84(0.51) 


-0.86(0.63) 


-0.049(0.139) 


0.028(0.154) 


10.29(4.12) 


11.16(0.76) 


Aq-F-5 


0.25 


0.39 


-1.05(0.59) 


-0.87(0.68) 


-0.002(0.131) 


-0.024(0.184) 


9.44(5.59) 


9.76(3.89) 


Aq-G-5 


0.59 


0.35 


-1.03(0.58) 


-1.02(0.93) 


0.019(0.126) 


0.043(0.177) 


11.05(1.89) 


10.47(2.13) 


Aq-H-5 


0.54 


0.46 


-0.79(0.49) 


-0.87(0.67) 


-0.056(0.175) 


0.148(0.420) 


10.88(2.10) 


11.82(0.88) 



Table 6. Chemical properties of the outer diffuse stellar haloes. The first column shows the encoding name of the simulations, Fp rog and 
F prom are the fractions of stars formed in the progenitor systems and of stars which have been part of a galactic wind, a and Age are 
the mean velocity dispersion (in km s" 1 ) and age (Gyr)of the stars. [Fe/H] and [O/Fe] are median abundances. Standard deviations are 
included within parentheses. The labels "prog" and "sat" denote the values corresponding to stars formed in situ and to stars acquired 
from satellites. 



Halo 


-Fprog 


■fprom 


[Fe/H]P r °s 


[Fe/H] sat 


[0/Fe]P r °s 


[0/Fc] sat 


AgeP r °s 


Age sat 


Aq-A-5 


0.35 


0.20 


-2.04(1.25) 


-1.44(1.31) 


0.018(0.168) 


0.145(0.203) 


10.79(3.93) 


11.08(2.14) 


Aq-B-5 


0.13 


0.12 


-1.59(1.39) 


-1.70(1.12) 


0.019(0.154) 


0.119(0.307) 


10.14(2.43) 


9.28(2.58) 


Aq-C-5 


0.19 


0.21 


-2.04(1.46) 


-1.31(1.14) 


0.054(0.149) 


0.089(0.205) 


11.54(0.92) 


11.29(1.11) 


Aq-D-5 


0.19 


0.24 


-1.76(1.24) 


-1.36(1.07) 


0.018(0.160) 


0.100(0.249) 


10.59(2.47) 


11.08(1.22) 


Aq-E-5 


0.25 


0.20 


-1.66(1.17) 


1.39(1.16) 


0.019(0.155) 


0.079(0.242) 


10.82(1.52) 


11.32(1.67) 


Aq-F-5 


0.03 


0.23 


-2.11(1.65) 


-1.32(1.09) 


0.019(0.160) 


0.019 (0.302) 


11.14(6.99) 


9.64(3.41) 


Aq-G-5 


0.25 


0.05 


-1.78(1.44) 


-1.85(1.22) 


0.018(0.161) 


0.089(0.320) 


10.55(6.45) 


10.54(4.07) 


Aq-H-5 


0.15 


0.29 


-1.73(1.64) 


-1.21(1.02) 


0.022(0.136) 


0.075(0.222) 


11.80(2.36) 


11.21(1.21) 



consequence of the bursty and extended star formation ac- 
tivity which gives rise to this component, allowing new stel- 
lar populations to form from material enriched by SNIa. Our 
SN feedback model allows us to track gas particles which 
have been promoted from the cold to the hot phase, thus 
becoming part of a galactic wind. We have checked how 
much of each disc formed from gas that was ever part of a 
wind, finding this to vary between ~ 10 and ~ 30 per cent. 
Hence, the great majority of disc stars formed from material 
that was acquired directly by infall and was enriched in situ 
without cycling through a galactic fountain. 

In central spheroids the level of a-enrichment increases 
with in situ fraction as indicated by the median abundances 
in Fig. [7] (lower panel). This is also consistent with the age 
distribution of these stars which, in general, are old, forming 
in strong starbursts which are too short for SNIa enrichment 
to be important (Fig. [3J. We find that more than 40 per 
cent of the stars in our central spheroids formed from gas 
that had been promoted at least once, showing that these 
dense and strongly bound regions are able to retain much 
of their SN ejecta. Since these regions are a-enhanced, the 
cooling and re-accretion of the promoted gas particles had 
to occur in less than » 1 Gyr. 

The outer and inner stellar haloes are more a-enhanced 
than the central spheroids. For these components, the frac- 
tion of stars formed from gas that had participated in a 
wind varies from 5 to 50 per cent, with the largest frac- 
tions of such stars in the inner haloes. We found that these 



'recycled-material' stars are primarily accreted rather than 
in situ objects, so were typically formed in satellite galaxies 
which were later disrupted. This is a consequence of our SN 
feedback model which results in larger promotion rates (rel- 
ative to the star formation rate) in systems with shallower 
potential wells: smaller haloes experience more significant 



winds (jScannapieco et al. 2008; Sawala et al.|20T0 De Rossi 
et al.|2010||Sawala et al.||2011[ ). 



So far we have analysed the abundance distribution of 
our stellar components without regard to where the stars 
formed. However, there are different levels of a-enrichment 
in the three spheroidal components, even though all are 
dominated by old populations with quite similar star for- 
mation histories. This suggests that the assembly history of 
each component plays an important role. Hereafter, we focus 
on this point. 

Properties of the in situ and accreted subpopulations of 
each component in each Aquarius galaxy are listed in Tables 
[3]to [6] Each table refers to a different component and gives 
for each galaxy the stellar fractions formed in situ and from 
previously promoted g well as the median [Fe/H] and 

[O/Fe] and the mean stellar age separated into in situ and 
accreted components. 

The in situ stars in discs have younger and more diverse 
ages than those of any other (sub)population, as can be seen 
from the penultimate column in Table [3] In particular, the 
few disc stars acquired from disrupted satellites are not only 
older but formed, on average, over shorter time intervals. As 
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Figure 1. Distribution of e = J z /J z , max (S) as a function of energy E for our eight Aquarius haloes. The last panel shows the color code 
used to represent the [Fe/H] abundances of star particles. Energies have been zero-pointed at the value for a circular orbit at the virial 
radius. The energies which separate our three spheroidal components are indicated for galaxy system in the relevant panel. For the sake 
of clarity, only 10% of the particles are displayed. 



a consequence, these stars tend to have lower metallicity and 
to be more a-enhanced than the in situ disc population. In 
situ disc populations also tend to have lower velocity disper- 
sions than accreted disc populations (see columns 4 and 5 of 
Table|3|. This suggests that accreted stars should contribute 
primarily to the metal-poor end of a thick disc component 
( Abadi et al.|2003| [Wyse|2010| . 

Our central spheroids are populated primarily by old 
stars which formed in situ over a relatively short time inter- 



val (~ 1.5 to 4 Gyr; see Table[4]|. The exceptions are Aq-E-5 
and Aq-F-5, which have been significantly rejuvenated by 
younger stars, as indicated by their larger age dispersions. 
Generally, accreted stars in the central spheroids tend to be 
older, to have formed over very short intervals, to have lower 
metallicities and to be more a-enhanced than those formed 
in situ. 

Our inner stellar haloes are formed by a mixture of in 
situ and accreted stars (Table [5J. In general, the latter have 
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Figure 2. Histograms of stellar mass as a function of e = J z /J z 
satellites have been excluded. 



t(E) for our Aquarius haloes within twice r ga i at z = 0. Stars in 
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Figure 3. Cumulative stellar mass as a function of age for stars in the disc (green), the central spheroid (black), the inner halo (violet) 
and the outer halo (red) for all simulated galaxies. 



similar metallicities but higher a-enhancements. They also 
tend to be older and to have formed over a shorter period 
of time. Similar trends hold for stars in the outer haloes 
(Table [6|, although in this case accreted stars tend to be 
more metal-rich as well as to have higher a-enhancements 
than the in situ population and there is no clear age sepa- 



ration between the two populations. Apparently, at least a 
small fraction of outer halo stars formed in situ from SNIa 
enriched material. These may also have formed at lower red- 
shift, giving rise to the large age dispersions listed for several 
galaxies in Table [6] Recall that stars formed in situ account 
for less than 20 per cent of the stellar mass in the outer 
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Figure 4. Distributions of [Fe/H] (upper raw) and [O/Fe] (lower raw) for stars in the discs (green), central spheroids (black), inner 
stellar haloes (violet) and outer stellar haloes (red) of four Aquarius galaxies. 
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Figure 5. Distributions of [Fe/H] for stars in the discs (green), central spheroids (black), inner stellar haloes (violet) and outer stellar 
haloes (red) of Aq-E sampled at three different levels of resolution, declining from Aq-E-5 to Aq-E-7. 



haloes and that all bound satellites have been removed and 
are considered as separate objects. Hence, we are here fo- 
cussing on the "smooth" component of stellar haloes. Our 
stellar haloes have a dual formation process since the ratio 
of accreted to m situ stars typically increases from 1 to 5 



between the inner and outer haloes (see also Zolotov et al. 
( 2009 1). We will explore this issue in more detail in the next 
section. 

Old metal-rich stars can be viewed as chemical fos- 
sils ( jBland-Hawthorn fc Freeman|2003| . We have measured 
the fraction of stars both older than 11 Gyr and with 
[Fe/H]> in each component of each of our simulated galax- 
ies. As expected the largest fractions are found in the cen- 
tral spheroids, ranging from eight down to one per cent. The 
other components have very few stars with these characteris- 
tics (fractions always below 1%). This is consistent with the 
fact that only in deep potential wells is it possible to reach 



high levels of enrichment in such a short period of time. Al- 
though we find no correlation with F plo& , we did detect a 
trend with the strength of the first starbursts: those cen- 
tral spheroids with a large fraction of old metal-rich stars 
have more than 40 per cent of their total stellar mass al- 
ready formed at z > 3 (i.e. older than 11 Gyr). Conversely, 
the other central spheroids were able to form less than 25 
per cent of their final stellar mass by this time. The low- 
est fraction of old metal-rich stars is found in Aq-B-5 which 
made only 15 per cent of its final stellar bulge by z ~ 3.5. 
Hence, the fraction of old and metal-rich stars in the cen- 
tral spheroids is linked to how fast star formation proceeded 
in the first massive structures and so is of interest for con- 
straining galaxy formation models, as claimed in previous 
work (e.g. Wyse|20"To and references therein). 
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Figure 6. Accreted (dashed) and in situ (solid) fractions as a function of formation redshift for stars in each of our four components 
and for all eight of the Aquarius galaxies. Each of the panels corresponds to a different component as labelled, while the colours of the 
curves identify the individual galaxies. The key is: Aq-A-5 (black). Aq-B-5 (red), Aq-C-5 (magenta), Aq-D-5 (cyan), Aq-E-5 (green), 
Aq-F-5 (blue), Aq-G-5 (violet) and Aq-H-5 (light green). 



4 THE INFLUENCE OF THE MASSES OF 
DISRUPTED SATELLITES ON THE 
ABUNDANCE PATTERNS OF ACCRETED 
COMPONENTS 

The abundance patterns of the various components of our 
simulated galaxies depend on the fraction of their stars 
which was accreted. However, the dispersions in these re- 
lations are large, suggesting that other factors are also im- 
portant. In particular, the masses of the disrupted satellites 
vary from halo to halo (e.g. Cooper et al.||2010|) and these 



differences could be reflected in the chemical patterns. In 
order to investigate this issue, we estimate the percentage 
(fmassivc) of the total stellar mass in a given dynamical com- 
ponent formed in satellites with stellar masses larger than a 
certain value. As can be seen in Fig. |8j the stellar mass of 
contributing satellites varies widely. The variation is partic- 
ularly marked for the inner and outer stellar haloes which 
were assembled primarily from satellite debris. As a refer- 
ence value for Fmassive we take that measured for stellar mass 



5 x W g MQh , similar to the stellar masses of the surviving 
satellites at z — in the Aquarius simulation^] 

In Fig. [9] we show the mass fraction of stars formed 
in satellites more massive than 5 x IO^q/i" 1 (F m 



massive J 



a function of F prog . There is a clear anti-correlation which 
simply reflects the fact that the larger the mass fraction 
formed in satellites, the larger the chance that massive ob- 
jects have contributed. In the discs and central spheroids, 
which formed primarily in situ, the accreted stars come 
mainly from small satellites (i.e. at most ~ 20 per cent from 
massive ones). The central spheroids of Aq-E-5 and Aq-F-5 
show the largest accreted fractions and a substantial fraction 
of their stars come from massive satellites, ~ 0.30 — 0.40 per 



3 Note that we do not use the ratio between the mass of the pro- 
genitor and that of the satellite because mctallicitics are related 
to the total stellar mass which could have enriched the system 
and to its potential well which affects the overall star formation 
efficiency and the impact of outflows. 
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cent, consistent with the fact that they alsos show the lowest 
median [O/Fe] (Fig[7|. The outer and the inner haloes have 
the largest contributions from stars accreted from massive 
satellites while the lowest such contributions are found in 
discs. As we show in Fig. [8] the mass functions of the con- 
tributing satellites vary significantly from galaxy to galaxy, 
leaving their imprints on the chemical abundances. 

In Fig. [ToJ we show the [Fe/H] and [O/Fe] medians as 
a function of F ma ssive. The stellar populations of the inner 
and outer haloes tend to be more metal-rich when larger 
fractions of their stars are accreted from massive satellites 
(see Fig. 10 upper panel). In this case, the correlation is 
clearer than that with the fraction of in situ stars. The inner 



haloes of Aq-B-5 and Aq-F-5 show similar levels of enrich- 
ment to their outer haloes. From Fig. [8] we see that the outer 
and inner haloes of these particular galaxies have similarly 
important contributions from stars formed in massive sys- 
tems, while for the other galaxies, the contribution of stars 
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Figure 8. Fraction by mass of each stellar component acquired 
from satellites more massive than a given stellar mass as a func- 
tion of this stellar mass. The dashed lines denote the reference 
stellar mass, 5 X lO 9 M0h _1 See Fig[6]for the colour code. 



acquired from massive satellites decreases significantly as 
one moves from the inner to the outer stellar halo. Thus, 
the stellar masses of the satellites where accreted stars were 
formed is an important factor in determining their mean 
level of enrichment: more metal-rich halo stars tend to come 
from more massive satellites. 
As shown in Fig. 



10 (lower panel), a-enhancement 



tends to increase with a decreasing contribution from mas- 
sive satellites for all three spheroidal components. The low- 
est median [O/Fe] values are found in central spheroids. 
Inner haloes extend over a larger range of both [O/Fe] 
and Fmassivc, while outer stellar haloes show the highest a- 
enhancements, reflecting the fact that they are dominated 
by old stars formed during short starbursts. 

A dual origin for halo stars is suggested by observa- 
tions of the Milky Way which show the inner halo to have 
a stellar metallicity distribution peaking at [Fe/H]~ —1.6 
while the outer halo is dominated by less enriched stars with 
peak [Fe/H]« -2.2 (|Carollo et ~aT1|2007[). There are also 



new observations suggesting differences in kinematics and 
a-enrichment for stars in the Milky Way stellar halo which 
might indicate the presence of two distinct stellar popula- 



tions (e.g. Nissen fc Schuster 2010 Beers et al. 20111. We 
now explore if the relative abundances of the inner and outer 
stellar field haloes of the Aquarius systems are comparable 
to these observed values and if a radial trend is a common 
characteristic of all our simulated systems. We stress the 
fact that in this section, we are analysing the diffuse stel- 
lar population in each component. Surviving satellites are 
explored in the last section of this paper. 

We have estimated the differences in median [Fe/H] be- 
tween the inner and outer stellar haloes of each of our eight 
simulated galaxies and looked for relations to the fraction of 
in situ stars. In all systems, the outer halo is less chemically 
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Figure 9. Fraction by mass of each stellar components accreted 
from satellites with stellar masses larger than 5 X W 9 MqIi~ 1 
) as a function of F prog . See Fig[7]for the symbol code. 



enriched than the inner one, with differences varying from 
r 0.6 dex to ~ -1.0 dex (Fig. 



111. The values quoted 



by |Carollo et aL] ( |2007[ ) for the Milky Way imply a differ- 
ence of ~ —0.6 dex, consistent with the range in the simu- 
lations where these gradients reflect the gradients in assem- 



bly history. Zolotov et al. (20091 and Zolotov et al. (20101 



reached similar conclusions, but used a different approach, 
so caution is needed in any detailed comparison with our re- 
sults. The dispersion in chemical gradients reflects a variety 
of factors, including the accreted mass fraction, the mass 
function of disrupted satellites, the ages of their stars, the 
mixing effiency, and the strength of outflows. For example, 
Aq-G-5 has the largest chemical differences between the in- 
ner and outer haloes due to a larger contribution of debris 
from massive satellites to the inner halo than to the outer 
one. a-enhancements depend principally on the star forma- 
tion timescales and the level of chemical mixing both in 
disrupted satellites and in the main progenitor systems. The 
inner haloes of our Aquarius galaxies are typically ~ 0.8 dex 
less cv-enhanced than their outer haloes, but in Aq-A-5 and 
Aq-C-5 the inner haloes are slightly more a-enriched than 
their outer components for given AFprog. Fig. [8] shows that 
these two inner haloes obtained their accreted stars from 
much smaller satellites than the corresponding outer haloes, 
contrasting with the other Aquarius systems. 

A similar comparison can be done for the central 
spheroids and the inner haloes. In this case we found metal- 
licity differences ranging from 0.9 to 0.4 dex with central 
spheroids always more enriched than the inner haloes. Based 
on observations reported in the literature (e.g WysepoTo I , 
we estimate a corresponding difference of ~ 0.35 in the Milky 
Way, in reasonable agreement with our results even though 
none of our simulated systems has a morphology resembling 
that of the Milky- Way. These trends are shown in Fig |ll| 
Larger differences in metallicity, are again associated with 
larger differences in formation history between the two com- 
ponents. The central spheroids are also systematically less 
a-enriched than the inner haloes, and again this trend is 
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more important for larger differences in the way the compo- 
nents were assembled. 



5 SURVIVING SATELLITES 

The chemical abundances of Milky Way satellites and nearby 
dwarf galaxies have been a subject of much research. Ob- 
servations suggest significant differences between the abun- 
dance patterns of stars in the Galactic halo and those in 
dwarf galaxies: stars in dwarfs tend to have lower [a/Fe] 
than similar metallicity stars in the Galactic halo (e.g. Tol- 
|stoy et al.|[2003t > even at low [Fe/H] levels. Thus it appears 
that the Galactic halo could not have been built from sys- 
tems similar to the current satellites. In this Section, we 
compare these two populations in our simulated systems. 
For this analysis, we have only considered satellites repre- 
sented by more than 2000 total particles. Two of our galaxies 
(Aq-E-5 and Aq-H-5) are excluded from this analysis, since 
they have no satellites satisfying this numerical requirement. 
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of the differences between the fractions of stellar mass formed in 
situ (AFp r0 g). See Fig. [7] for the symbol code. 



In Fig. |12[ we show cumulative stellar mass as a func- 
tion of stellar age for each satellite of the remaining six 
galaxies. Clearly, star formation histories differ considerably, 
with most showing several starburst episodes. For compar- 
ison, we have included the star formation histories of the 
outer haloes of the corresponding systems. Most surviving 
satellites have experienced several starbursts and are domi- 
nated by younger stellar populations than the corresponding 
outer haloes. Few satellites have significant old stellar popu- 
lations and only one halo, Aq-B-5, has satellites dominated 
by older stars than those in the corresponding outer halo. 
This figure shows the variety in the star formation histories 
of satellites inhabiting similar haloes, illustrating the diffi- 
culty of comparing the Milky Way to any specific simulation. 
For this reason, we compare the chemical properties of ac- 
creted stars in our four stellar components to those of stars 
in the surviving satellites on a statistical basis. 

Fig. " 



median [Fe/H] for stars in surviving satellites (solid lines) 
and for accreted stars in the outer haloes (dotted lines). 
Accreted stars tend to be a-enhanced, reflecting the fact 
that they formed quickly and in short bursts, so that SNII 
enriched them but SNIa exploded too late to affect their 
abundances. In contrast, surviving satellites contain signif- 
icantly more weakly a-enhanced stars at all metallicities. 
This is a consequence of their extended star formation his- 
tories (Fig. 12 I. As shown in previous papers, our feed- 



back model successfully drives mass-loaded galactic winds 
out of star-forming galaxies with a wide range of masses, 
and with a specific mass-loss rate which depends on the po- 
tential well depth of the systems. These winds regulate star 
formation activity and determine the fraction of available 
baryons which are turned into s tars (jScannapieco et al.| 2008[ 



20091 jj2010| |Sawala et al.|2010||DeTiossi et al.|2010||Sawalal 
et al.|| 2011[ ). As a result, the chemical properties and ages 



of the stars in surviving satellites are statistically different 
in our simulations from those of stars acquired by satellite 
disruption. 

An interesting point of comparison between our cur- 
rent model and observation concerns the percentage of very 
metal-poor stars ([Fe/H] < —3) in the surviving satellites. In 
our simulations this varies from 0.5 to 20 per cent (with an 
average of ~ 10 per cent), depending on details of individ- 
ual star formation histories. These stars are a-enhanced and 
old. The apparent lack of such low-metallicity stars in real 
systems has been suggested as a possible problem for hier- 



archical galaxy formation models in the past (e.g. Shetrone 
et al.|2001||Helmi et al."|20 06). However, more recent studies 



13 shows the distributions of median [O/Fe] and 



have discovered extremely metal-poor stars in dwarf galax- 
ies ( |Kirby et al.|2008[|2009"l|Norris et al.|2008||2010[ ). These 
new observations, together with improved metallicity indi- 
cators (ptarkenburg et al. 2010), are causing a revision of 
ideas about the low metallicity tail of the stellar populations, 
which are clearly more extensive than previously thought. 
These tails are of great interest for investigations of the 
possible impact of other processes such as pre-enrichment 
by Population III stars which might have led to a floor for 
Population II metallicities, leading to a deficit of extremely 
metal-poor stars in dwarfs. Improved observational surveys 
of low metallicity stars in nearby satellites are thus very 
much relevant for constraining galaxy formation models. 

Our results are encouraging in that they show differ- 
ences between the chemical properties of accreted halo stars 
and those of stars in surviving satellites which are quite 
similar to the trends observed in the Milky Way. Higher res- 
olution simulations are clearly needed to confirm these find- 
ings, but they already suggest that these trends are not only 
compatible with the standard model of hierarchical forma- 
tion in a ACDM cosmology, but may, indeed, provide direct 
evidence for the assembly paths it predicts. 



6 CONCLUSIONS 

In this paper we study the chemical patterns of the stellar 
populations in eight Milky- Way mass-sized galaxies simu- 
lated within the ACDM scenario. Our version of GADGET- 
3 includes the modelization of SN feedback describing the 
release of energy and chemical elements in the interstellar 
medium so that the chemical enrichment of baryons can 
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Figure 12. Cumulative stellar mass as a function of age for surviving satellites (violet dashed lines) and for the outer stellar haloes 
(black solid lines) of six of our Aquarius galaxies. Aq-E-5 and Aq-H-5 have only very low-mass satellites and have not been included in 
this plot. 



be followed consistently as galaxies are assembled. This is 
the largest sample of high resolution Milky- Way mass-sized 
galaxies so far run with a physical motivated SN model fact 
that provides us with the unique possibility to distinguish 
those chemical properties which are common to all galaxies 
from those that are determined by the particular history of 
assembly of each galaxy. 

Our results can be summarized as follows 

• Discs: The surviving discs show very similar chemical 
properties. They are more chemically enriched than the stel- 
lar halo but slightly less enriched than the central spheroids. 
The stars in discs show the lowest levels of a-enhancement. 
This is consistent with being formed from SNIa pre-enriched 
material. These stars formed mainly in situ. Most of the 
stars in the discs formed from material which has never been 
part of a galactic wind. In general, the discs formed inside- 
out as shown by |Scannapieco et al.| (2009, 2011[), and during 
this process different starbursts enriched continuously the 
gas. We identified that up to 15 per cent of the stars in the 
discs formed in satellite systems confirming previous results 
( Abadi et al.|2003 1 . These stars tend to older and formed in 
shorter starbursts in comparison to the dominating stellar 
population formed in situ. Accreted stars have systemati- 
cally higher a-enhancement and are more supported by ve- 
locity dispersion. These stars could contribute to a thick disc 
component. Most of the accreted stars come from satellites 
smaller than 5 x 10 9 Mmfo -1 total stellar mass. 



• Central spheroids: They are formed mainly by old stars, 
born in starbursts with different duration ranging from 2 
to 4 Gyrs approximately. The stellar populations in central 
spheroids have been formed in situ with contributions of less 
than ~ 20 per cent of stars formed in satellites. Two systems 
out of eight (Aq-E-5 and Aq-F-5), show larger contributions 
of accreted stars. These central spheroids have been also 
recently rejuvenated by news stars formed mainly in situ 
by the infall of enriched gas (mean stellar age of the new 
population 4-5 Gyrs). This can be confirmed by the lower 
a-enhancement showed by these two central spheroids with 
respect to the rest of the systems. Very old and metal-rich 
stars are principally located in the central spheroids and 
their frequency can be correlated with rate of star formation 
in the earliest epoch of formation. 

• The inner and outer field haloes: As one moves out- 
ward the contribution from disrupted stars to the field stellar 
haloes increases dramatically. The inner haloes have contri- 
bution from both kind of stars, but the outer haloes have 
been mainly formed by accreted stars. The mean metallic- 
ity shows a large dispersion which correlates with the frac- 
tion of stars acquired from massive satellites (i.e. larger than 
5 x 1O 9 M0/i _1 ) so that the larger this fraction, the larger 
the mean level of enrichment of the halo. This is valid for 
both the inner and the outer haloes. 

We find that in all simulated stellar haloes the stars in the 
outer region are systematically less enriched. The differences 
in the chemical properties between the inner and the outer 
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Figure 13. Distributions of median [O/Fe] (upper panel) and of 
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haloes (dotted lines) and for stars in surviving satellites (solid 
lines) . 



field haloes are not only the consequence of the different 
relative contributions of in situ and accreted stars ( |Zolotov] 
et al. 2009 1 but also reflect the mass function of the satellites 



which contributed with stars to form them. 

The level of a-enhancement is higher in outer field stel- 
lar haloes reflecting the fact that they have been mainly 
assembled from old stars formed in smaller structures, in 
single short-duration starbursts. Only stars accreted from 
very massive systems tend to be, on average, less a-enriched 
at a given metallicity than the global population of the halo 



iZolotov et al. 20101. The masses of the specific satellites 



where the accreted stars were formed are reflected in the 
chemical properties of the inner and outer haloes, imprint- 
ing features which distinguish one system from another. 

• Surviving satellites: Each simulated Milky-Way mass 
galaxy has a system of surviving satellites with particular 
characteristics. The surviving satellites show a variety of 
star formation histories but in general, they formed stars 
in a series of starburst episodes which extended even to the 



present time. As a consequence, their chemical patterns dif- 
fer from those of the stars that contributed to the formation 
of the stellar components of the main galaxies. On average, 
stars in surviving satellites are younger, low metallicity and 
low a-enriched. 

Although there are many open questions and issues to 
be solved regarding galaxy formation, particularly those that 
could be related to the formation and surviving of disc galax- 
ies, it is encouraging that the global chemical trends ob- 
served in nearby galaxies can be obtained in ACDM scenar- 
ios. Our work shows the importance of studying chemical 
evolution to understand galaxy formation and the relevance 
of doing so within a cosmological context. 
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